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Abstract 
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1 Introduction 



In this article we continue the investigation of the finite-size scaling properties of reaction- 
diffusion models defined on an one-dimensional lattice of L sites. We consider the example 
of a special model with only one kind of particles denoted by A in which the elementary 
processes are diffusion, coagulation and decoagulation. This model is integrable and allows 
the occupation probability <rii> (t) of site i to be computed analytically A detailed analysis 
for periodic boundary conditions has been given in our previous article |J (hereafter referred to 
as I) , where we focused our attention on the finite-size scaling properties of the concentration: 

c(L,t) = }E <ni>{t) (1.1) 
L i=i 

in the scaling limit L — > oo, t — > oo keeping the scaling variable z = fixed. We computed 
the corresponding scaling functions exactly and observed that they do not depend on the 
initial conditions. 

In this paper we consider the same model with open boundary conditions. In particular 
we want to understand how the boundary conditions modify the finite-size scaling properties 
of the concentration c(L,t). 

Let us first describe the model and outline the method which allows us to compute the 
concentration. The model is defined by three elementary processes: 

• diffusion rate D + A <-> A + 

• coagulation rate c A + A — > A (1.2) 

• decoagulation rate a A —>■ A + A 

where A and denote an occupied and an empty site, respectively. The empty interval 
approach, which we will use to calculate the concentration, requires that coagulation and 
diffusion take place with equal rates. Therefore we have to choose c = D. Furthermore it 
is convenient to fix the time scale so that D = 1. The system is now controlled only by 
one parameter, namely the decoagulation rate a. Since all rates are chosen to be left-right 
symmetric the model is invariant under space reflection. In paper I we summarized how 
the master equation of the reaction-diffusion system defined by ( |L2| ) can be mapped onto a 
one-dimensional quantum chain 0. In this formulation the time evolution of the model is 
described by the euclidean Schrodinger equation: 

||P) = -H\P) (1.3) 

where the Hamiltonian H is given by: 

H = J2Hi (1.4) 

i=i 

and: 

Hi = + 2E1 1 ® Egi + (1 + a) {E™ ® + E} 1 <g> E™ x ) (1.5) 
- ® + E}° ® E^x + ® EH, + ^ ® E^) 
-aiEl^El^+El^E™,). 
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Here are 2x2 matrices denned by (E° ' )k,i — &k&i acting on the i th site. Notice that 
open boundary conditions are already implemented in Eq. ( |1.4| ). Our aim is to compute the 
mean value of the concentration: 

c(L,t) = j £ <m> (t) = j E(0|^ n exp(-^)|P Q ) (1.6) 
L i=i L i=i 

for an arbitrary initial state \Pq) (for notations we refer to paper I). 

In our first paper we have shown that for a ^ the Hamiltonian H with periodic boundary 
conditions can be mapped onto the well-known anisotropic Xl^-chain with periodic boundary 
conditions. Since the corresponding similarity transformation is local, it also applies to the 
case of open boundary conditions. One obtains the Hamiltonian: 

h xy = ~E (Vx+i + rVM +1 + < + *Ui - v-^ 1 ) (i.r) 
z i=i 

where: 

TJ = y/\ + a (1.8) 

is a real parameter. Notice that although the original system ( |1.4|) is defined with open 
boundary conditions, in Eq. ( |1.7| ) we do not obtain a quantum chain with free boundaries since 
there are additional fields -(erf + o~ z L ) at the ends of the chain. The XF-chain is integrable, i.e. 
it can be diagonalized exactly in terms of free fermions, which allows to write the Hamiltonian 
(|1.7|) in the diagonal form [|J: 

L-l 

H XY = A fc 4a fc , (1.9) 

fc=0 

where a\ and at are fermionic operators: 

{4, a,} = 5 ktl . (1.10) 
The Afc are the corresponding fermionic excitation energies and are given by: 

Ao = 0, (1.11) 

(irk \ 
T) + r]- 1 - 2 cos — J , fc = l,...,L-l. (1.12) 



The excitation with vanishing energy (|1.11|) for arbitrary t] is related to a hidden quantum 



group symmetry M and implies the levels of the spectrum to be at least two-fold degenerated. 
The second equation ( 1.12j) is just the dispersion relation of the Xl^-chain and we observe that 



the system is massless for 77 = 1 (a = 0). Taking all combinations of the excitation energies 
Afc into account, one can construct the spectrum of H XY and therewith the spectrum of H. 

One way to compute the concentration c(L, t) would be to diagonalize H directly. This 
turns out to be very complicated since it is not known in advance which of the eigenvectors 
enter the expression for the concentration and thus all of them have to be calculated. Therefore 
we prefer a different approach based on so-called empty intervals probabilities instead of spin 
configurations || |(| (for a compact summary we refer to paper I, the most general formulation 
of this approach can be found in [0). The time evolution of the probability n, t) of finding 
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an empty interval of length n extending from site j ' + 1 — ^ to site j + | is described by a 
closed system of linear differential equations. Since 1 — fl(j — ^, 1, t) is just the probability to 
find the j th site occupied the concentration c(L, t) is simply given by: 



c(L,t) = yX;<n,->(f) = ^(l-n(i-J,M)). (1.13) 

A detailed derivation of this formalism can be found in paper I. Here we present a modified 
version of the differential equations for open boundary conditions. If the interval does not 
touch the boundaries we have: 



dfl(j,n,t) 



m V 2 (p(j - ~, n + 1, t) + Q(j + ~ n + 1, tj) (1.14) 

+ - + Q(j + ^,n-l,t)-2{ V 2 + i) n&M) > 

where | < j < L — | and n = 1, • • • , L — 1. If on the other hand the interval touches the 
boundaries one obtains the 'boundary equations': 

^ = ry 2 Q(— ,n+l,t) (1.15) 

+ nC-g-, n - 1, t) - (v 2 + 1) n, t) , 

dQ(L-%,n,t) 9 . n + l , 

1 m 2 ' = V 2 n(L-—,n + l,t) (1.16) 

Th — 1 Th 

+ n(L — , n - 1, t) - (rj 1 + 1) tt(L - -, n, t) , 



dQ(±L,t) 



2 

0. (1.17) 



9* 

As in paper I in all equations we use the convention: 

n(j,o,t) = i. (Lis) 

Therefore this system of equations is inhomogeneous. A particular solution is obviously given 
by: 

n{j,n,t) = 1 (1.19) 

corresponding to an empty lattice. In this notation the space reflection is described by the 
map: 

P : (j,n)»(L-j,n). (1.20) 

As one can easily check the system of differential equations (|1.14 )- (|1.17|) is invariant under P. 
Notice that, in opposition to the case of periodic boundary conditions where the Fourier 
transformed system has only L — 1 degrees of freedom, this system has L(L + l)/2 degrees 
of freedom. It grows quadratically with lattice length L and thus it is more difficult to solve. 
The central mathematical aspect of this paper is the exact solution of this problem. 

Because of the reduced number of degrees of freedom compared to the original physical 
system the set of differential equations ( 1.14Q - (|1.17[ ) can only cover a part of the full dynamics 
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of the Hamiltonian (|1.4j ). As we will see in Sec. ^| this part corresponds just to those states 
of the XY chain with only one or two excited fermions. The advantage of the empty interval 
formalism is that c(L,t) can be computed easily via Eq. (|1.13|) . 

This article is organized as follows. In Sec. |^ we solve the system of differential equations 
(PID -( P7|) . We then derive exact expressions for c(L, t) for special initial conditions in 
Sec. |3| and analyse the corresponding finite-size scaling properties. Finally we summarize our 
results. 



2 Solution of the System of Differential Equations 



In this section we solve the system of differential equations ( |L14|) -( |I.17|) . For this purpose 



we first determine the complete set of solutions <j)A(j,n,t) of the homogeneous eigenvalue 
problem: 

d 

—<p A (j,n,t) = -A(f) A (j,n,t) , (2.1 
where Eq. ( |1.18| ) has to be replaced by: 



<MJ, 0,t)=0 (2.2 

in order to impose homogeneous boundary conditions. It will be convenient to introduce the 
notations: 

n n 

where x and y take the values < y < x < L. In this notation the space reflection reads: 

P : {x,y)^{L-y,L-x). (2.4 
For < y < x < L the system of differential equations ( |1.14| )-( |i.l7| ) now has the form: 

d(p(x,y,t) 



dt 



d<f)(x, 0, t) 
dt 

d(f>{L,y,t) 
dt 

d<f>{L,0,t) 
dt 

For 0<x = y<Lwe have: 
or: 



r l 2 ( ( p(x + l,y,t) + ( p(x,y-l,t)) (2.5 
+ (j>{x - 1, y, t) + <j>{x, y + l,t)- 2{if + 1) <j>{x, y, t) , 

ri 2 <P(x + 1, 0, t) + <j)(x - 1, 0, t) - (r] 2 + 1) 0(x, 0, t) , (2.6 

t/ 2 0(£, V-l,t)+ <t>{L, y + ht)- ( V 2 + 1) <j>(L, y, t) , (2.7 

. (2.8 

<f){x,x,t) = , (2.9 



<f>(x,x,t) = l (2.10 

if we choose the system to be homogeneous or inhomogeneous, respectively. For easy reference 
we illustrated the organization of the complete system for L = 6 sites in Fig. 1. Here the 
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L(L + l)/2 boxes represent the degrees of freedom <f>(x,y). The number in each box denotes 
the strength of the diagonal contribution in the differential equations, where d — — (2 + a) — 
— (1 + 1] 2 ), while the arrows denote the flow of probability to the neighbours with respect to 
the couplings 

l + a = i] 2 . 

The zeros on the hypotenuse x = y represent the contributions from the inhomogenous part 
of the system. Since we consider only the homogenous equation due to Eq. ( [2.9|) these terms 
have to be neglegted throughout this section. 

Before calculating the solutions of the eigenvalue problem of Eq. (|2.1|) we give a short 
description how to proceed. Looking at Fig. 1 we find that the boxes on the short sides of 
the triangle form a decoupled (2L — l)-dimensional subsystem of differential equations. As a 
consequence we expect two kinds of solutions: 

• The first kind is obtained by solving the subsystem separately. Since each solution 
imposes inhomogeneous Dirichlet boundary conditions it can be extended uniquely to 
the interior of the triangle. This procedure yields a set of 2L — 1 solutions of Eq. (|2.1| ) 
forming a complete set of eigenfunctions for the short sides of the triangle. 

• The second kind of solutions is now calculated using the completeness of the first kind 
solutions. They are forced to vanish on the boundary and therefore to obey homogeneous 
Dirichlet boundary conditions. 

Let us now study the solution of Q2.1|) in more detail. For the moment we choose rj = 1 (see 
Eq. (|1.8| )). Then the differential equations acting on the short sides of the triangle are: 



X 

D 2 y ML,y) = -AML,y), (2.12) 



D*<f> A (x,0) = -A<j> A (x,0), (2.11 

l2 



where x, y = 1, • • • , L — 1. Here D 2 and D 2 are discrete second derivative operators in x and 
y direction, respectively. D 2 for example is defined by D 2 (f)(x, y) = (f)(x + l,y) + 4>(x — l,y) — 
2(j)(x,y). 

The probability for an empty lattice 4>a(L, 0) decouples completely and remains constant 
in time. Thus 4>\(L, 0) takes nonzero values only if A = 0. In this case the second derivative 
of 0a=o vanishes and 0o becomes a linear function of x and y, namely <f>o(x,y) ~ x — y due 



to Eq. (|2.9|) . For A ^ the condition (j) A (L,0) = gives an additional constraint which 



determines the solutions of Eqs. (|2.11 ) and ( |2.12j ) to be simple oscillations 



, , . . nkx . . . irky . 1 

(p(x, 0) ~ sin — - , (p(L, y) ~ sm — - , (2.13) 

( 7r ^ 
Au = 2 l — cos — 

V L 

where k — l, ... ,L — l. The values of A give the spectrum of eigenvalues. Comparing it 
with the spectrum of the XF-chain ( |l.9| ) we recognize that it corresponds to one-fermion 
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excitations. For a given 7^ every linear combination of the solutions of Eq. ( |2.13| ) solves 
Eqs. ( 2.11 ) and (|2.12 ). Thus these eigenvalues are twofold degenerated. However, we can 



choose the functions <pk to be eigenfunctions of the space reflection P (see ( |2.4|) , too. It will 
be convenient to introduce the notation P(f>^ = ±(— l) fc 0^. So far we have determined the 
required 2L — 1 solutions on the short sides of the triangle which still have to be extended to 
the interior. Here the differential equations act according to (D^ + -D^)0 = — A0 where the 
operator on the left side is just the discrete Laplacian. The solutions are given below in Eqs. 
( P7T5D and (|2~T6| ). 



Now we turn to the calculation of the remaining ( L ~ 1 K L ~ 2 ) solutions of the interior of the 
triangle. These solutions are plane waves ( [2.17| ) vanishing on the boundary of the triangle. 
The corresponding eigenvalues A = + A; are identified as two-fermion excitations of the 
Xy-chain. 

Summing up all results we have three types of solutions: 



• The stationary solution: 

One-fermion excitations (k = 1, . . . , L — 1): 



Mx,y) = ^ T 1 - (2.14) 



, , , . 1 / . irk . nk \ 

= -j= ^ sm ~y x - sm -j-y) , ( 2 - 15 ) 

1 / 2y nk 2x irk \ 



• Two-fermion excitations (k, I = 1, . . . , L — 1; k < I): 

2 / nk Til ttI irk \ 

® kl {x,y) = — I sm —x sm —y — sm —x sm —y J . (2.17) 

This is the complete set of L(L + l)/2 solutions of the homogeneous system. For rj ^ 1 we 
obtain deformations of these solutions. In order to avoid too much details, we only want to 
present our results: 



The stationary solution: 



Mx,y) = ; ' _ 2L (2-18) 



1 — rj~ 



with eigenvalue A = and parity P = +1. 
One-fermion excitations (k = 1, . . . , L — 1): 



4(x,y) = ^ T ^^±^)sin^a;-(^±^)sin^) (2.19) 
with eigenvalue A& = 77(77 + r]~ l — 2 cos ^) (see Eq. ( |1.12| )) and parity P = ±(— l) k . 
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• Two-fermion excitations (k, I = 1, . . . , L — 1; k < I): 

2 / irk ttI ttI irk \ 

9ki{x,y) = —r] y I sin — x sin — y — sin — x sin — y I (2.20) 

with eigenvalues A k i = r](2r] + 2ri~ l — 2 cos ^ — 2 cos ^) and parity P = (—l) k+l+1 . 

It is easy to check that these solutions reduce to Eqs. ( |2.14| )-( p.l7| ) in the limit rj — > 1. 

Let us now study the orthogonality relations of the different types of solutions. The 
functions (j) k , fc ^ form a complete system of eigenfunctions on the short sides of the 
triangle except the point (L,0). The corresponding scalar product denoted by ( | )boundar y 
takes only these points into account. For arbitrary functions f,g we define: 

(/ | ^boundary = E ^/C*. 0) + £ T, 2 ^ f {L, y)g(L, V ) . (2.21) 

The solutions (p. 19 ) are normalized in such a way that: 

(0fc I 0fc')boundary = (2.22) 

where fc, fc' = 1, . . . , L — 1 and a, (3 = ±. For the functions (|2.20 ) we define a scalar 
product in the interior of the triangle: 

(/ 1 s>buik = E E v 2(x - y) f(x, y)g{x, y) ■ (2.23) 

x=l y=l 

Then the functions <&ki (|2.20|) obey the relation: 

(®ki I ®k'i')buik = hk'Sw ■ (2.24) 

These relations are independent of the value of rj. Unfortunately the solutions cj)^ take nonzero 
values in the interior of the triangle. Therefore the scalar product {<&ki | 0fc/)buik is not equal 
to zero. These projections have to be taken into account if one expands an arbitrary function 
Q(x, y) in terms of the eigenfunctions as we will do in the next section. 



3 Finite-Size Scaling of the Concentration 

In this section we use the complete set of eigenfunctions found in Sec. |2| in order to 
calculate the time evolution of the concentration for an initially fully occupied lattice and 
an arbitrary value of r\. In the case of rj — 1 the finite-size scaling limit is computed and 
compared with the corresponding result for periodic boundary conditions. 

Before calculating the concentration for special initial conditions we consider the general 
case. Let VLq{x, y) be an arbitrary initial choice of empty interval probabilities satisfying 
condition (|1.18|) . The time evolution then takes the form: 

n(x,y,t) = l + u M^y) + J2T,^rk(x,y)e- A '' t + J2 E u; k ^ kl (x,y)e-^ +A ^ (3.1) 

k=l a=± k=l l=k+l 
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where the constant 1-function is just the particular solution ( |1.19| ) and Afc is defined in Eq. 
( |1.12| ). The coefficients to , and tu^i are obtained by expanding il (x,y) in terms of the 
eigenfunctions (|2.18| )- (|2.2CI| ). Therefore we first consider Eq. (|3.1|) for t = at the point (L, 0). 
Since all eigenfunctions except </>o vanish at this point we have: 



wo = n (£, o) - 1 



(3.2) 



Then we consider the short sides of the triangle, where the eigenfunctions $ki vanish. Here 
the coefficients can be obtained using the relation ( |2.22|) : 



± 



J k ~ \0fc I ^0 — 1 — ^0 00 ) boundary • 

Finally we calculate the remaining coefficients Uki with the help of Eq. Q2.24| ): 



L-l 



U k l = (®M | ^0 - 1 - ^O0o)bulk "EE ^k'i^kl | 0fc')bulk • 

k'=l a=± 



(3.3) 



(3.4) 



We now apply this general procedure to the case of an uncorrelated initial state where the 
occupation probability p is equal for all sites (see paper I). Then we have: 

n (x, y ) = (i- P y-y . (3.5) 

Straightforward but rather laborious calculations show that the concentration is given by: 



c(t, L) 



1 - n~ 2 

hrV-(i-P) L ) 



1 — 77 

+ L^'- ^ £ (1 " (1 " P) V(-l) fc )(^ + ~ 2(-l) fc ) 



(3.6) 



V 



k=l 



X 



sin 2 (^) / I 



A* 



+ 



L-l L-l gin 2 rrk gin 2 ^ 



L 3 r] 



E E 

k=l l=k+l 
k+l odd 



cos ^ — cos j; 



A fc A fc 
1 



cos ^ — COS y 



exp (-A fe t) 



1 1 

Afc A, 



1 1 

Afc Ai 











(Afc 


Afc/ 







AfcA/ H A^Afc 



where Afc is defined by: 



A* 



1] (1 — p)r/ + (1 — p) rj —2 cos 



exp (-(Afc + Aj)*) 



T 



(3.7) 



Now we consider the finite-size scaling limit of the concentration. As pointed out in paper I we 
assume that in the finite-size scaling regime L — > oo, * — > oo, z = 4| fixed the concentration 
behaves like: 

c(z, L) = L x [F (z) + L~ y F(z) + ■•■]. (3.8) 



S 



In this formula Fq(z) denotes the scaling function and x the scaling exponent. L y F(z) is 
the leading correction term. The functions F (z) and F(z) are supposed to depend only on 
z. Since the finite-size scaling hypothesis is only valid for critical systems corresponding to 
massless quantum chains we have to restrict ourselves to the case rj — 1. In this case the 
concentration takes the form: 



c(t,L) 



1 - (1 -p) L 
L 

4 sin 2 2£ sin 2 f r 

E — 1 



(3.9) 



fc=0 1=1 
k+l odd 



COS 



7rfc 



COS 



■kI 
L L 



COS ^ - COS ^ 



1 



1111 

Afc A; A fc A; 



1 



^a-aiJU-x;.)' 1 -^- 1 '' 



2 2 
^ + 



exp(-(A fc + A/)i) 



A fc A, A,A fc _ 

Performing now the finite-size scaling limit we get for the functions F (z) and F(z): 



oo oo 



w = i+Se e 



fc 2 + z 2 



7T" 



/, u 1=1 O' 
k+l odd 



cxp 



Z7T 



-(fc 2 + / 2 )) 



(3.10) 



oo oo 



F(z) 



fe=0 J=l 
fe+i odd 



-7T 



, (tf + /2 )(A; 4 + /4) g (1 + 6 i_ ? )( A; 4 + /4) + m 



(k 2 ~l 2 ) 2 



(k 2 - l 2 f 



(3.11) 



x exp 



ZTt 



(A; 2 + f 



and for the exponents x, y: 



x = — 1, 



y = 2 



(3.12) 



The functions F and F can be written in a simpler form by rewriting the sums in Eqs. (|3.10|) 
and ( |3.11| ) in terms of // = and z/ = ^y^: 



F (z) = 1 - 7^0(2) 5 

7z z 1 / / N 2 

^0) = -^^1(^)^2(2) - — g G {z)gz{z) - — ( <7o(z)<fe(z) - (01(2 



-(1 + 6 



1 — p. /I 



-fl-oO^O) + 2(^1^ 



(3.13) 



(3.14) 



where the functions gi(z) are defined by: 



9i{z) 



E 



„- S o**0x+i)' eXP 



Z7T 



2 ( ^ + 2 ; 



9* 



3. 



(3.15) 
(3.16) 
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Especially we have the identity: 

0i(*) = 2 (o,^) (3.17) 

where 8 2 (u,t) is a Jacobi-Theta function. Now we are able to compare the finite-size scaling 
behaviour for open and periodic boundary conditions. In both cases we find the scaling 
exponent x — — 1 and the correction exponent y = 2. Thus the critical exponents are not 
influenced by the boundary conditions. The scaling function and the correction function on 
the other hand are different as can be seen in Fig. 2 and Fig. 3. However we find that in both 
cases the scaling functions and the correction functions are related to Jacobi-Theta functions. 
As in the case of periodic boundary conditions we observe that the scaling function ( |3.10| ) is 
independent of the initial probability p whereas the correction function ( p. 11 ) is not. 



4 Conclusions 

Based on the empty interval probability approach we solved the coagulation-decoagulation 
model defined on a lattice of L sites with open boundary conditions. In contrast to ordinary 
mean-field calculations our results take the microscopic fluctuations of the system into ac- 
count. We computed the concentration c(L,t) exactly and investigated its finite-size scaling 
properties. Our aim was to compare the results with those of our previous paper 0, where 
periodic boundary conditions have been considered. 

Mathematically the time evolution of the empty interval probabilities is described by a 
system of L(L + l)/2 linear differential equations (in opposition to only L — 1 equations in 
the case of periodic boundary conditions). It is the first time that such a system has been 
solved in the theory of reaction-diffusion processes, and we sketched the general method how 
to find the solutions. We obtained three types of solutions wich can be interpreted as states 
with none, one and two fermions in an excited state. 

We used the expression for c(L, t) to compute the finite-size scaling expansion of the 
concentration and compared the result with the case of periodic boundary conditions. It 
turned out that one obtains the same type of expansion (with the same powers in 1/L). 
As in the periodic case we found that the scaling function is independend of uncorrelated 
homogeneous initial conditions. However, they are different for open and periodic boundary 
conditions. On the other hand the correction function is influenced by the initial condition 
similar to the periodic case. In this paper we only investigated homogeneous uncorrelated 
initial conditions. Initial states with small clusters are studied numerically. The results are 
presented in ||. It turns out that the scaling function remains unchanged even for that kind 
of initial conditions. 

Using the empty interval approach we can compute the local concentration Cj(t). As 
mentioned in the introduction, this formalism does not cover the full dynamical space of the 
physical system. In particular it is impossible to compute correlations between concentrations 
located at distinct points. For this reason we want to outline a possible generalization. An 
empty interval can be understood as a two-point object. Therefore it is natural to consider 
also generalized 2n-point objects consisting of n separated empty intervals. These objects are 
also expected to be described by a set of linear differential equations as well. In particular 
they are nothing but n-point correlation functions of the occupation number if we choose the 
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lengths of all intervals to be equal to one (similar to Eq. ( |1.13| )). Of course we expect a higher 
number of fermionic excitations to be involved. 
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Figure 1: Structure of the system of differential equations for L = 6 sites. The figure is 
explained in the text. 
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